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Abstract. 

The Brownian motion of a harmonically bound quantum particle and coupled to a harmonic 
quantum bath is exactly solvable. At low enough temperatures the stationary state is non-Gibbsian 
due to an entanglement with the bath. This happens when a cloud of bath modes around the particle 
is formed. Equilibrium thermodynamics for particle plus bath together, does not imply standard 
thermodynamics for the particle itself at low T. Various formulations of the second law are then 
invalid. First, the Clausius inequality can be violated. Second, when the width of the confining 
potential is suddenly changed, there occurs a relaxation to equilibrium during which the rate of 
entropy production is partly negative. Third, for non-adiabatic changes of system parameters the 
rate of energy dissipation can be negative, and, out of equilibrium, cyclic processes are possible 
which extract work from the bath. Conditions are put forward under which perpetuum mobile of the 
second kind, having several work extraction cycles, enter the realm of condensed matter physics. 



Introduction. There are not two fundamental theories of nature, quantum mechanics 
and thermodynamics. There is only one: quantum mechanics, while thermodynamics 
must emerge from it. The universal character of equilibrium thermodynamics led to the 
general expectation that in one way or another, thermodynamics will apply to the full 
quantum domain [jl]]. Few people have taken the painful road to check this emergence, 
yet this is what we have set out to do. Here we discuss the results for quantum Brownian 
motion that have been presented [Q, ^ and were discussed in the scientific literature [Hp. 

Brownian motion has numerous applications in condensed matter physics [^J, |[ |7|, 8], 
atomic physics [Q], quantum optics and chemistry [£p. Some realizations involve weak 
coupling with the thermal bath [|9|] . However, there are well-known experimental situa- 
tions, which are essentially far from the weak-coupling regime. Here standard thermo- 
dynamics may not apply. The main example of this is the case of weak links between 
superconductive regions, the so-called Josephson junctions, in their overdamped regime 
[ pTj| , P"Tj |, where the relevant ranges of parameters were achieved already twenty years 
ago. Even in quantum optics, which has often been satisfactorily described by weak- 
coupling theories [Q], there are recent experiments showing the necessity for moderate 
and strong coupling approaches [fl"2|]. 

The Hamiltonian. We consider an 'ideal gas' of non-interacting harmonic oscillators 
coupled to a bath. For the total Hamiltonian Jt? tot = J$? + Jt? B + M'j we thus assume ^\ 
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with ffl describing the particle, <%e B the bath, and the interaction. 

The bath is assumed to have uniformly spaced modes CQ ( - — iA, / = 1,2,3,---. The 
thermodynamic limit for the bath is taken by sending A — > 0, which induces relaxational 
behavior. For the couplings we choose the quasi-Ohmic Drude-Ullersma spectrum, 
where [§] J(co) = j7r^(c?/m ( .«o.)5(c()-fi).) = ya>T 2 / '(co 2 +T 2 ), where yis the damping 
constant. We shall assume that the Debye cutoff T is large. 

We describe an ensemble of closed total systems with conserved energy, except for 
the periods when work is done on it by externally changing m or a. We take as initial 
density matrix the Gibbs distribution exp(— /3^ ot ) /Z tot , with Z tot the partition sum. 

Equilibrium state. In Gibbsian equilibrium of the total system the free energy reads 
F tot (r, 7) = F g (T, 7=0)+ F p (a, 7, T, m, T) . The free energy of the bath itself, 

T l"^ 7T 2 T 2 

Fb(Ti7=0) = -1 (tota (l-.-f) = -_, (2) 
is of order 1/A, showing the extensivity of the bath. The Brownian particle adds to this 
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It contains three characteristic, temperature independent of frequencies. For small 7 
(underdamping) they read: co 1 2 = ±ico Q + J /2m, Gt) 3 = T — y/m. On the other hand, 

for strong damping y 2 3> am, co l = a/y, co 2 = (y/m) (1 — am/ y 2 ), ct) 3 = T — y/m. 

Effective temperatures. We shall now study two objects, T x = a(x 2 ) and T p = (p 2 )/m, 
that would in classical equilibrium be equal to T and which we shall interpret below 
as effective temperatures. As in the classical situation, it holds that T x = 2adF p /da, 
T p = —2mdFp/ dm. For 7^0 one gets the weak-coupling result known from textbooks, 

U = \T X + jT p = T X = T p = ^d) coth^/3/?(0 , with co Q = \fajm. At large T one gets 
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At low T and for strong damping one has 



The fact that T p > T x > at T = is related to the quantum nature of the problem. 

Generalized thermodynamic formulation. The Wigner function has a quasi- 
Gibbsian expression, since there occur two temperature-like variables, 

-p 2 /2mT p -ax 2 /2T x 

W(p,x) = W p (p)W x (x) = - 7 == r—— . (6) 

y/2nmT p ■ s /2%T x /a 

The Boltzmann entropy of momentum and coordinate is 

f , x r / x 1—^ 1 emT„ 1 eT x 

S p = -j dpW p (p)ln[W p (p)Vh] = -\n-^, Sx =-lnj±, (7) 



while the von Neumann entropy reads [||] S vN = (v+ 1) ln(v + 5) — (v — j) ln(v — 5), with 
v = ApAx/ h = \J (p 2 ) (x 2 ) /h = ^mT p T x / Ti^fa. The first terms in its large v-expansion, 
S vN = lnv + 1, coincide with the total Boltzmann entropy S B = S p + S x . 

Internal energy and interaction energy. The energy of the central particle reads U = 
(Jj?) — jT p + jT x . The energy of the cloud of bath modes that surround the particle is 
U m = U m - U B (y =0)-U = U p -U = TdF p /dT. At high temperatures one has U = 
r + /3/z 2 (2a + yr)/24m, U int = ptfyT/lAm. Since the energy of the cloud involves h, it 
is a quantum effect. At low temperatures one has, if T is large, U- t = hy/2nm + <ff(T A ). 

Implementation of the first law. Given the Hamiltonian J4?(p,x) of the subsystem, 
its energy is U = (Jf?) = J dpdxJ4?(p,x)W(p,x), where W is the Wigner function 
of the subsystem. When we change a system parameter, the energy changes as dU = 
d / J^W = fJ^dW + fW dJt? = d^ + dW. The first term represents the variation due 
to the statistical redistribution of the phase space, which we identify with the change in 
heat d£}. The last term results from the change in the Hamiltonian, so it is a mechanical, 
non-statistical object, which we associate with the work dW done by external sources. 

Generalized free energy and the second law. The definition of the effective tempera- 
tures admits a thermodynamical interpretation. The free energy F for a two-temperature 
system is defined as [ |13| ] F = U — T p S p — T X S X = — jT p \nmT p — jT x \n(T x /a). For 
adiabatic changes in m or a one has dF = —S x dT x — S p dT p + dW , with work dW = 
(dJ4?/dm)dm+ (dJff /da)da = —T p dm/2m + T x da/2a. Due to the first law this yields 
the second law for situations with two temperatures, d£2 = T p dS p + T x dS x , in close anal- 
ogy with those proposed recently for nonequilibrium glassy systems [ [T3| ] and black 
holes [JH]]. Notice that F pertains to the particle alone, and does not satisfy Gibbsian 
thermodynamics, while the Gibbsian F p of Eq. @, relates to the whole equilibrium sys- 
tem, i.e., to the particle and the cloud of bath modes around it. There are many physical 
systems, such as a Josephson junction strongly coupled to the electromagnetic field, 
where the natural object to study is F, relating to properties of the junction only. 

Violation of the Clausius inequality at low T. The Clausius inequality dJ2 < TdS vN 
is one of the formulations of the second law. At T = it says that no heat can be taken 
from the bath, at best heat can go from the subsystem to the bath. 

When we change a — > a + da, d£} is of order — T 2 da, while TdS vN is of order — Tda. 
In the case da > 0, where an amount of work dW rew ~ da is done on the system, the 
Clausius relation is thus violated at low T. Likewise one can consider the variation of 
the (effective) mass m. Here one has d£} = hydm/lnm 2 . Now there is a transfer of heat 
even if the bath temperature is zero. Thus, violation of the Clausius inequality is even 
stronger in this case. This situation with d W KW < corresponds to the work performed 
by the system on the environment. That the heat comes from the cloud of bath modes, is 
confirmed by the fact that d^ rev = — dC/ int . 

Violation of the Landauer bound for information erasure. A further aspect is the 
squeezing of phase space and entropy, relevant for computing in the quantum regime. We 
have shown that the so-called Landauer bound — di2 > kT \nl for the minimal energy 
dispersion when erasing one bit of information is violated in a similar manner [£p. This 
connection arises because the Landauer bound is just based on the Clausius inequality. 

Dynamics from a non-equilibrium state. Energy oscillation at low T. We consider 
the dynamical evolution of a system initially in equilibrium characterized by a spring 
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FIGURE 1. a) Normalized excess energy f { as a function of at/y, for different values of the temper- 
ature. Upper curve: T — ► °°. Middle curve: r = Tia/2y. Lower curve: at T = it is non-monotonous, 
b) Schematic plot of the cyclic changes in the spring constant, where successive cycles are slower 
and slower, h characterizes the size of the change and t denotes the dimensionless time. The interval 
— °° < t < t j marks the process that establishes the nonequilibrium state at t = 1 1 . The picture shows three 
full cycles, in the intervals t ( <t < t j+l , (i = 1,2,3). Their start and end points are indicated by bullets. 



constant a , which at t = is instantaneously changed to a { = a. These parameters are 
connected as a = (1 — (X Q )a and we assume that \ a Q \ <C 1. For strong damping one has 
for the energy of the subsystem U(t) = U(°°) + (haa /2ny)f l (at /y). The behavior of 
fl(o) at different temperatures is presented in Fig. la. For a Q > it says that, after 
initially energy has been put on the particle by the change of a Q — > a > a Q , this energy 
leaks away into the bath. At low T, however, too much leaks away, and a part has to 
come back. This "bouncing" is familiar of the noise correlator, which behaves like f x . 

Entropy production. The rate of production of Boltzmann entropy is found after a 
lengthy identification of its flux. The result can be negative at moderate T [Q]. For T = 
and strong damping one has with o = at jy and o ■ = 0.879: djS x /dt ~ +(Xq(o — 
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+ e), where e = am/y ■ This is negative for a n 



e < a < o„ 
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m/y, the characteristic timescale of the momentum. 
Smooth changes of the spring constant. Let us now consider the case where, starting 
from the equilibrium state a(— °°) = a, the spring constant a{t) = [1 — a(t))a, is slightly 
changed (|ce(f)| 1) in a smooth way. The rate of work added to the system is dW /dt — 
dW rev /dt + dTl/dt, with the adiabatic ('reversible') and dispersed contributions 
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where X = yt/2m and C x a function like /j of Fig. la. Integrating over the full change 
(the whole region where a' ^ 0) one confirms that the dispersion for a completed, cyclic 
change (oc- = a.) of system parameters is nonnegative, the Thomson formulation of the 
second law [JT5J]. It is also nonnegative for noncyclic (a- ^ a^) but completed changes. 

Energy dispersion. If <x(x) = (X m h(Q.t), with £1 a slow rate of change, and 7 is large, 
expansion in £1 is possible. Omitting numerical factors one gets the structure 



dn 

dT 



hai&ti 



yT 

ha 



2r>2 



y L £l 



(9) 



The full integral is always positive and for moderate T the first term is dominant. How- 
ever, for small T the term h'h" may imply that dfl/d? is negative. This is a firm statement, 
also valid when starting in equilibrium, since work is the energy added to the total sys- 
tem. It is the more surprising since the relevant domain the characteristic timescale of the 
change of a, I/O. < h 2 / (jT 2 ), may exceed all other timescales, including the quantum 
timescale h/T, so that applicability of thermodynamics was to be expected. 

Perpetuum mobile with many work extraction cycles. The h'h" term can cause extrac- 
tion of work. Let there be JV non-overlapping Gaussians, h(x) = exp( — \x 2 ), as depicted 
in Figure lb. Cycles are the intervals t t <t < t i+l (i > 1). Each new cycle is slower than 
the previous one, £l n , j < Q. n . One can make a total number JV ~ 1/T of cycles with 
equal yield, where jV is parameterized by v, and total yield 

ha rPK&i/Zn dyy ^ _ TZO^yT 2 rP^/to V y 

IkjT Jpm r jin v+I 1 y + I 3 y 3 ' tot 3 ha Jpno. v /2n v + I x y + I 3 y 3 ' 

(10) 

with I { = \-Jk, 7 3 = 3 y/H/A. The minus sign of W tot indicates that work is performed by 
the system on the environment. This is possible because, in order to make the extraction 
cycles, one had to start from the equilibrium state a(— 00 ) = and change a up to a{t x ). 
In this first part the energy dispersion was IT(? 1 ) = +(X^hyQ\/2ATZa > \W tot \. 

Feasibility. The harmonic oscillator can be interpreted as an LC circuit [jn5p. x may 
correspond to the charge Q on a capacitor, 1 /a to its capacitance C, m to an inductance 
L, p to a phase 4>, 7 to a resistance R, and r\ (t) to a random electro-motoric force. In this 
setup there should be nothing difficult in varying L ~ m or C ~ I /a. We have proposed 
to test the violation of the Clausius inequality in mesocopic linear circuits [jlTp. This 
amounts to measuring (8Q 2 ), (54> 2 ) and the produced work. The conditions to do this 
experiment were reached 20 years ago, and one of the quantities, (SQ 2 ), was measured 
already and agrees perfectly with the theoretical predictions [|T8], [19]]. 

Has thermodynamics been violated? Let us recall that our results hold also for Af 
non-interacting Brownian particles in a bath. Our conclusion is that thermodynamics 
does not always work when, in the quantum regime, one considers the Brownian particle 
in its reduced Hilbert space, thus summing out the bath. We should admit that at low 
enough T the interaction energy between system and bath becomes relevant, since the 
damping constant is fixed. In a very strict definition of thermodynamics one may claim 
that there is no reason why thermodynamics had to apply. 

For this reason our surprising findings of the breakdown of the Clausius inequality can 
in principle be viewed as results outside the domain of applicability of thermodynamics. 
However, the negativity of the Boltzmann entropy production sets in already at moderate 
T, more or less in the same domain where energy initially put on the particle starts to 
dissipate to the bath in a non-monotonic fashion, even in limit of large damping. 

Our second effect in that regime, the presence of many work extraction cycles ("per- 
petuum mobile") involves only the energy budget of the total system, and should be 
more surprising. However, when starting in an equilibrium state and counting all the 



energy, the total dispersion is positive, as demanded by an exact theorem Q15|]. 

These extraction cycles relied on the fact that the rate of energy dispersion by the total 
system can be negative. This is a firm and unexpected statement about the total system, 
that may have started in equilibrium. Thus positivity of the rate of energy dispersion 



[ |20| ] is not always a good formulation of the second law. 

The new aspects of quantum Brownian motion arise from quantum entanglement: A 
complete description in terms of a wave function is possible only for a closed system; 
subsystems are necessarily in a mixed state. Thus the quantum Gibbs distribution is 
not an adequate candidate for the description of the quantum subsystem non-weakly 
interacting with its thermal bath. Connections of these findings with NMR physics [|21|], 
mesoscopic work sources Jill ] and with Maxwell's demon [ p2| ] are discussed elsewhere. 
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